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HIGH-TEMPERATURE THERMAL STORAGE SYSTEMS 
FOR ADVANCED SOLAR RECEIVERS 
MATERIALS SELECTION 

D. F. Wilson, J. H. DeVan, and M. Howell 
ABSTRACT 

Advanced space power systems that use solar energy and Brayton or 
Stirling heat engines require thermal energy storage (TES) systems to 
operate continuously through periods of shade. The receiver storage units, 
key elements in both Brayton and Stirling systems, are designed to use the 
latent heat of fusion of phase-change materials (PCMs). The power systems 
under current consideration for near-future National Aeronautics and Space 
Administration space missions require working fluid temperatures in the 1 1 00 
to 1400 K range. The PCMs under current investigation that have liquidus 
temperatures within this range are the fluoride family of salts. However, 
these salts have low thermal conductivity, which causes large temperature 
gradients in the storage systems. Improvements can be obtained, however, 
with the use of thermal conductivity enhancements or metallic PCMs. In fact, 
if suitable containment materials can be found, the use of metallic PCMs 
would virtually eliminate the orbit associated temperature variations in TES 
systems. 

The high thermal conductivity and generally low volume change on 
melting of germanium and alloys based on silicon make them attractive for 
storage of thermal energy in space power systems. An approach to solving 
the containment problem, involving both chemical and physical compatibility, 
preparation of NiSi/NiSi 2 , and initial results for containment of germanium 
and NiSi/NiSi 2 , are presented. 


1. INTRODUCTION 


Thermal energy storage (TES) systems are required for continuous operation of 
advanced space power systems that use solar energy and Brayton or Stirling heat 
engines. Storage systems have been designed to use the latent heat of fusion of rnetals 
or salts. Energy is stored during the heating and melting of a phase-change material 
(PCM) in the sunlit portion of the orbit. The PCM cools and freezes, thereby releasing 
heat to the engine, during the eclipse period of the orbit. In a latent-heat system, heat 
is theoretically absorbed and liberated at a constant temperature. However, in actuality, 
a constant temperature is difficult to achieve because the solid phase of the PCM (lower 
thermal conductivity) forms on the heat-transfer surfaces. 1 The current reference PCM 
thermal energy storage system for the Brayton Cycle is based on LiF-CaF 2 eutectic 
contained in Haynes-1 88 canisters. Characteristics of this system are listed in Table 1 . 
Possible salt PCMs for the Stirling engines are listed in Table 2. 
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Table 1 . Reference storage system for 
Brayton Cycle Solar Dynamic systems 


PCM 


Material 

LiF 20 CaF 2 

Melt temperature 

1040 K 

Heat of fusion 

Specific 

81 5 kJ/kg 

Volumetric 

1.71 kJ/cm 3 

Thermal conductivity 

Solid 

4.0 W/nrK 

Liquid 

1,6 W/rn-K 

Density, liquid 

2.1 g/cm 3 

Void formation 

22% AV/V, 

Container 

Material 

Haynes-1 88 

Maximum temperature 

~1100 K 

System 

Specific heat storage 9 (maximum) 

300 kJ/kg 

Volumetric heat storage 9 (maximum) 

1 .3 kJ/cm 3 


9 Estimated using the dimensions of a heat storage 
canister proposed for a space station system, assuming 
complete melting and freezing of the PCM in each 
cycle. 

Source: Wichner, R.P.; Schmidt, M.A.; and 

DeVan, J.H., Oak Ridge National Laboratory, "Improved 
Performance Options for Solar Dynamic Heat Storage 
Systems," personal communication to NASA Lewis 
Research Center, Cleveland, Ohio, 1988. 
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Table 2. Salt PCMs for Stirling engines 



T m 

Ah 


(K) (kJ/kg) 

T in oDtimum ranae (1 100 K-1200 K) 


— m ■ 

NaF-24 MgF 2 

1 1 03 r *' 655 

LiF 

1120 

1087 

NaF-27 CaF 2 -36 MgF 2 

1178 

-520 

T above ODtimum ranae (1200 K-1300 K) 


MgF 2 -43 CaF 2 

1227 


NaF-36 MgF 2 

1260 


NaF 

1268 

789 

NaMgFg (line compound) 

1303 

711 

T below ODtimum ranae (1000 K-1 100 K) 


LiF-20 CaF 2 

1040 

815 

NaF-32 CaF 2 

1083 

-540 


Source: Wichner, R.P.; Schmidt, M.A.; and 
DeVan, J.H., Oak Ridge National Laboratory, 
"Improved Performance Options for Solar 
Dynamic Heat Storage Systems," personal 
communication to NASA Lewis Research Center, 
Cleveland, Ohio, 1988. 


Several parameters are important in the choice of a PCM. 2 The melting 
temperature of the PCM must be higher than the operating temperature of the chosen 
thermodynamic heat engine. The heat of fusion and density determine the mass and 
volume of PCM required to store a given amount of energy. The volume change on 
fusion (void formation) and thermal conductivity affect heat transfer and, thus, container 
design. Solomon 3 has demonstrated the dependence of the overall performance of the 
TES system on latent heat and thermal conductivity. He reported that doubling the latent 
heat of fusion of the PCM not only reduces the mass of the PCM required for energy 
storage by ~ 50%, but also reduces the percent variation (charging/discharging) in the 
energy stored. Also, he showed that doubling the thermal conductivity reduces the large 
oscillations in temperature of the outlet gas temperature by ~30%. 

In general, salts and their eutectics have low thermal conductivities. However, it 
is desired that only a small temperature drop occur across the PCM during thermal 
discharge to minimize variations in the temperature of the outlet gas into the power 
system. Therefore, large heat-exchange surface areas are required to minimize the 
temperature gradients in the energy storage and removal system. In addition, these salts 
undergo volume changes (as much as 40%) 1 during their phase changes, and voids are 
generated. These voids will likely have an adverse effect on heat transfer into and out 
of the TES container. Also, high stresses may occur if the expanding liquid phase 
cannot communicate with the voids. 
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Several means for improving the performance of thermal energy systems have 
been considered by Wichner et al.* These include 

1. higher heat storage temperature, thus increasing the thermodynamic cycle 
efficiency; 

2. higher latent heat, thus reducing the total mass of the receiver; 

3. higher thermal conductance, thus allowing more of the PCM to melt in a given time 
with a given temperature gradient; and 

4. reducing void volume, thus increasing heat transfer and reducing thermal 
stresses. 

All of these approaches can be implemented using metallic PCMs. 

Wichner et al. considered (1) congruently melting line compounds or elements in 
which a single phase is in equilibrium with a single liquid phase at its melting point; 
(2) eutectics, in which two solid phases are in equilibrium with a single liquid phase at 
the liquidus temperature; and (3) solid solutions, in which a single solid phase melts 
incongruently between the solidus and liquidus temperature. These types of materials 
are regarded as ideal because they in theory would release energy isothermally on 
solidification. Table 3 lists the materials suggested by Wichner et al. as candidate PCMs. 
The materials are divided into four groups, and within each group the members are 
expected to behavior similarly. However, many of the parameters of interest such as 
latent heat of fusion, thermal conductivity and expansion, volume change on 
transformation of state, density, liquidus temperature, and amount of undercooling are 
not known for the silicon-based alloys. 

Low-density metals are of particular interest as possible advanced heat-storage 
media because of their high thermal conductivity and generally low volume change on 
melting as compared to salts. However, the candidate metals are generally reactive with 
commonly used containment materials; as a result, the range of suitable materials is 
narrowed to ceramics with low connected porosity and small pore sizes. The 
compatibility problem for these TES systems was approached from a hierarchy that 
considered chemical reactivity; solubility; and physical compatibility, especially thermal 
expansion mismatch. 


2. MATERIALS 


2.1 PHASE-CHANGE MATERIALS 

Based on the analysis of Wichner et al., three candidate PCM materials were 
chosen for the initial screening work. These PCMs were germanium, a NiSi-NiSi 2 eutectic 
alloy (Group VIIIA - Si systems), and a Mg 2 Si-Si eutectic alloy (Group IIA - Si systems). 


*Wichner, R.P.; Schmidt, M.A.; and DeVan, J.H., Oak Ridge National Laboratory, 
"Improved Performance Options for Solar Dynamic Heat Storage Systems," personal 
communications to NASA Lewis Research Center, Cleveland, Ohio, 1988. 
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Table 3. Silicon-based alloys for thermal energy storage 




T m 

(K) 

AH m 

(kJ/kg) 


Group VIIIA 

- Si systems 


PdSi-Si 

Eutectic 

1143 

480 

PdSi 

Compound 

1363 

520 

NiSi-NiSi 2 

Eutectic 

1239 

640 

NiSi 

Compound 

1265 

750 

Fe 2 Si 5 -Si 

Eutectic 

1480 

1000 

FeSi-Fe 2 Si 5 

Eutectic 

1485 

900 

Fe 2 Si 5 

Compound 

1493 

1180 


Group VII • 

- Si systems 


Mn^-Si^ 

Compound 

1418 

1130 

ReSi 2 -Si 

Eutectic 

1398 

948 


Group IIA • 

■ Si systems 


Mg 2 Si-Si 

Eutectic 

1219 

800-1000 

Mg 2 Si 

Compound 

1358 

1100 

CaSi 2 -Si 

Eutectic 

1253 

1100 

SrSi 2 -Si 

Eutectic 

1273 

870 

SrSi-SrSi 2 

Eutectic 

1317 

660 

SrSi 2 

Compound 

1423 

850 

BaSi 2 

Compound 

1453 

630 


Ge-Si Solution systems 


Ge 

Element 

1211 

510 

Ge-0.2Si 

Solution 

1270-1410 

650 

Ge-0.4Si 

Solution 

1340-1510 

830 

Ge-0.6Si 

Solution 

1430-1580 

1130 

Ge-0.8Si 

Solution 

1560-1640 

1330 

Si 

Element 

1687 

1800 
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Germanium is commercially available. As will be discussed later, the NiSi-NiSi 2 eutectic 
can be readily made, but preparation of Mg z Si-Si eutectic has presented many 
unresolved problems. 

2.2 CONTAINMENT MATERIALS 

Thermochemical analyses were performed to estimate the chemical stability of 
possible containment materials with molten silicon alloys and germanium at 1100 to 
1400 K. The analyses utilized FACT (Facility for the Analysis of Chemical 
Thermodynamics, McGill University) and identified alumina, graphite, silicon carbide, and 
boron nitride as potential containment materials at the respective melting points of the 
referenced PCMs. 

Porosity of the containment materials was considered. Low porosity and small 
pore size are necessary to limit the ingress of the molten PCMs into the walls. Such 
ingress can lead to mechanical fatigue of the containment materials from repeated 
liquid/solid transformation. Also, segregation of alloy materials and, hence, loss of 
performance are possible. Low porosity graphite, chemically vapor deposited (CVD) 
silicon carbide, and CVD boron nitride meet these criteria. In addition, to minimize the 
open porosity at the surface that contacts the molten PCMs, the surfaces of some 
graphite crucibles were coated with pyrolyticaily deposited carbon (PyC). Wetting of the 
container by the molten PCMs was also evaluated because wetting will aid the ingress 
of the molten PCMs into the walls of the container. 

Finally, physical compatibility is also an important consideration. The containment 
materials must have sufficient mechanical strength or coefficient* of thermal expansion 
to mitigate the effect of the volume changes of the PCMs from phase transformation and 
thermal expansion. The thermal expansion coefficients 4 of a few materials are presented 
in Fig. 1 . It should be noted that the expansion of graphite varies over a wide range and 
can, at least in theory, be chosen approximately equal to that of germanium. Silicon 
carbide and boron nitride have expansions less than that of germanium. Thus, cooling 
of germanium to room temperature in these materials can be accommodated. However, 
heating to the operating temperature with a solid, snug-fitting mass of germanium could 
result in failure of the container. On the other hand, alumina has thermal expansion 
greater than that of germanium. Thus, cooling of germanium in this material should 
result in failure of the alumina container. However, because of the greater mechanical 
strength of alumina vs germanium, failure did not occur. 

Alumina was selected as the capsule material for compatibility testing of the silicon 
alloys (except as noted in the following), while both graphite and alumina were used to 
contain germanium. Although it was predicted that graphite would react with the silicon 
alloys to form silicon carbide, compatibility tests were performed nevertheless to evaluate 
whether the resultant silicon carbide reaction product would act as a barrier to extensive 
attack. Alumina, graphite, pyrolytic carbon on graphite, and CVD silicon carbide and 
boron nitride were obtained from commercial sources. 
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Fig. 1. Thermal expansion coefficients of germanium and selected 
containment materials. 


3. EXPERIMENTAL PROCEDURES 


3.1 PREPARATION OF NiSi-NiSL, EUTECTIC ALLOY 

The appropriate amounts of nickel and silicon for 38 wt % silicon were arc-melted 
and homogenized under an argon atmosphere. The molten alloy was drop-cast into a 
water-cooled copper mold under an argon atmosphere. Samples from three different 
positions along the length of the casting were evaluated using mass spectrometry, X-ray 
diffraction, and metallographic techniques. 

3.2 COMPATIBILITYTESTING 

Specimens of the candidate thermal storage material and the potential container 
material were placed in an "inert" alumina crucible (with cap). Alternatively, specimens 
of the thermal storage materials were placed in a crucible (with cap) of the potential 
container material. The alumina or container material crucibles were, as shown in Fig. 2, 
encapsulated under vacuum (1 O' 1 mPa) in quartz. The assembly was heated to 50 K 
above the melting point of the candidate storage material and maintained at this 
temperature for 1 68 h. Subsequent to exposure, the specimens were subjected to visual, 
metallographic, and scanning electron microscopic analyses. 





• N Si/M Si 2 with BN 
in Alumina Oudble 

•»** * • T \ 




Fig. 2. Test assembly consisting of an alumina crucible and 
cap encapsulated in quartz. 
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4. RESULTS AND DISCUSSION 


4.1 CHARACTERIZATION OF NiSi-NiSL, EUTECTIC ALLOY 


Chemical analysis for nickel and silicon in the NiSi-NiSi 2 alloy was performed by 
mass spectrometry, and the normalized concentrations are presented in Table 4. The 
designations SD1, SD2, and SD3 refer to the top, middle, and bottom of the cast 
specimen, respectively. The targeted composition was 38 wt % silicon. As is shown, 
there was some variation in composition through the length of the casting, and the 
average silicon content was higher than the targeted concentration. However, the 
composition was well within the NiSi-NiSi 2 two-phase field and had a melting point quite 
close to that of the eutectic. 

Optical micrographs of specimen SD4 from the mid-section of the casting are 
shown in Figs. 3 and 4. The specimen was metallographically prepared and etched with 
50% HCI/10% Gly/10% HN0 3 . A typical eutectic structure is evident with a slight excess 
of the NiSi 2 phase. Comparing Fig. 3 to Fig. 4, it is seen that the distribution of this 
excess phase is not uniform and, in places, quite small. The cracks suggest severs 
stress generation during solidification and cooling. 

The results of the X-ray diffraction analysis confirm the presence of both the NiSi 2 
and NiSi phases. Those phases are also indicated by the phase diagram. 5 However, 
certain features of the X-ray diffraction pattern warrant further discussion. A segment of 
the diffraction pattern produced using copper /c-alpha radiation of wavelength 1 .54059 A 
is shown in Fig. 5. The diffraction peaks corresponding to NiSi 2 were well defined. The 
average value of the lattice parameter calculated from the identified lattice planes is 
5.396 ± 0.008 A. In fact, the (333) and (440) reflections gave a lattice parameter of 
5.400 A. That is to be compared to 5.406 A for NiSi 2 (48.90 wt % silicon) 5 with a 
CaF 2 (Cl) type structure. The remaining diffraction peaks were quite broad. Thiu 
broadening renders the analysis more difficult and ambiguous. In this case, the lattice 
parameter values of a = 5.62 A, b = 5.18 A, c = 3.34 A for orthorhombic NiSi (32.37 
wt % silicon) 5 were used to calculate the diffraction angles. Comparisons were made 
with the broad peak positions. There was good agreement between the centroids of the 
broad peaks and the calculated diffraction angles. The broadening could be a result of 
stresses generated during casting of this material, a third phase, or a range of 
stoichiometry for NiSi. However, such a range is not indicated in the phase diagram, 5 
and stresses may have only affected the NiSi diffraction peaks. Thus, the presence of 
a third phase could not be absolutely ruled out. 


Table 4. Spectrographic analysis 
of NiSi-NiSi 2 eutectic alloy 


Element Specimen 


(wt %) 


SD1 SD2 SD3 


Si 

Ni 


40.9 44.5 42.3 

59.1 55.5 57.7 


Y210938 
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Fig. 3. Optical micrograph of NiSi-NiSL, specimen SD4, etched with 50% HCI/10% Gly/10% HNO 


Y210937 
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Fig. 4. Optical micrograph of NiSi NibL, specimen SD4, etched with 50% HCl/10% Gly/^0% HNO 
As compared to Fig. 3, there is a significantly smaller quantity of the excess pnase NiSi ? . 
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Fig. 5. Segment of X-ray diffraction pattern of NiSi-NiS^ eutectic, using 
copper /(-alpha radiation of wavelength 1 .54059 A. 








13 


4.2 COMPATIBILITY STUDIES WITH NiSi-NiSL 2 EUTECTIC ALLOY 

Exposures of NiSi-NiSi 2 eutectic alloy to pyrolytic carbon on graphite, alumina, and 
CVD silicon carbide and CVD boron nitride in alumina crucibles were performed at 
1 289 K. In all cases, on completion of the experiments, it was found that the crucibles 
had broken. However, the NiSi-NiSi 2 ingot was totally retained within the crucible and 
appeared to have begun solidifying with the crucible still intact (see Fig. 2). Neither the 
coefficient of thermal expansion nor the volume change (expansion) on freezing of 
NiSi-NiSi 2 is known. Therefore, it can only be speculated whether failure occurred at the 
liquid-to-solid transformation or after further cooling. Initial metallographic results of 
these short-term isothermal screening tests of pyrolytic carbon, CVD silicon carbide, and 
CVD boron nitride appear favorable from a chemical compatibility perspective. 

An optical micrograph of NiSi-NiSi 2 tested at 1 289 K for 1 68 h in a crucible of 
pyrolytic carbon of thickness 0.005 cm (0.0020 in.) on graphite is shown in Fig. 6. Note 
that there is a very thin layer of the pyrolytic carbon attached to the alloy with little 
evidence of reaction for this isothermal short-term experiment. An SEM micrograph of 
a portion of the pyrolytic carbon on graphite near a fracture site of the crucible is shown 
in Fig. 7. This surface was in contact with molten NiSi-NiSi 2 at temperature. Note there 
are indications that the pyrolytic carbon has been pulled away from part of the graphite 
surface. Energy-dispersive X-ray (EDX) analysis indicated that there were more nickel 
and silicon on the pyrolytic carbon surface than on the graphite surface. However, this 
analysis was complicated by the presence of "foreign" elements such as alumina, iron, 
calcium, and silicon in the graphite. 

As shown in Fig. 8, a negative meniscus was observed with CVD silicon carbide. 
This indicates a nonwetting system. An SEM EDX examination of the NiSi-NiSi^SiC 
interface indicated that one small area of the CVD silicon carbide was wetted by the 
NiSi-NiSi 2 eutectic alloy at the test temperature. This specimen was sectioned, and the 
area of interest was prepared for metallographic examination. Currently, nothing unique 
has been observed in that area. Also, the EDX analysis indicated that a quite small 
quantity of alumina was on the surrounding surface of the NiSi-NiSi 2 alloy. The quantity 
appears to result from a dusting of alumina. 

In Fig. 9, strips of CVD boron nitride contained within NiSi-NiSi 2 and germanium 
after compatibility testing are shown. Note that the NiSi-NiSi 2 eutectic alloy has shattered 
the alumina crucible (and is in one piece), while the germanium did not. A metallurgical 
examination of the NiSi-NiSi^CVD boron nitride interface indicated good chemical 
compatibility and no migration of the PCM into the boron nitride. An optical micrograph 
of the section is shown in Fig. 10. 

4.3 COMPATIBILITY STUDIES WrTH GERMANIUM 

The short-term exposures (168 h at 1283 K) of germanium to graphite, PyC, CVD 
silicon carbide, and CVD boron nitride all indicated good compatibility. There were no 
detectible chemical reactions and no migration of the germanium into the test 
containment materials. There was adequate physical compatibility with graphite and 
alumina as all the capsules remained intact on cooling to room temperature (see for 
example, Fig. 9). However, as Lauf 6 has demonstrated, depending on the grade of 
graphite, there can be seepage of germanium through the walls during thermal cycling. 
Seepage occurred with Stackpole graphite but not with Poco graphite. Further, modeling 
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Fig. 7. SEM micrograph of a portion of the pyrolytic carbon on graphite near the 
fracture site of the crucible. This surface was in contact with molten NiSi NiSi, 
at temperature. 
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Fig. 9. Optical photograph showing CVD boron nitride contained within NiSi NiSL, and 
germanium. Note that one crucible is shattered. 
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work on thermal performance of heat receiver concepts using germanium as the metallic 
PCM 7 has shown virtual elimination of orbit-associated temperature variations in the 
storage systems. 

4.4 Mg 2 Si-Si EUTECTIC ALLOY 

The appropriate amounts of Mg 2 Si and Si for 57 to 58 wt % silicon combined in an 
alumina crucible (with cap). The crucible, in turn, was encapsulated under 
vacuum (10' 1 mPa) in quartz. The assembly was heated to 1423 K and allowed to 
furnace cool. After cooldown the quartz envelope was found shattered, There we re no 
indications of alloying. There are questions as to the usefulness of this eutectic alley for 
use in thermal energy storage systems because of the high vapor pressure associated 
with Mg 2 Si. 


5. SUMMARY 


The study of the compatibility of containment materials with metal-based PCM 
involved considering chemical reactivity; solubility; and physical compatibility, especially 
thermal expansion mismatch. 

Thermochemical analyses were performed to estimate the chemical stability of 
possible containment materials with molten silicon alloys and germanium at 1100 to 
1 400 K. The analyses utilized FACT and identified alumina, graphite, silicon carbide, and 
boron nitride as potential containment materials at the respective melting points of the 
referenced PCMs. 

Graphite was confirmed experimentally under short-term isothermal conditions to 
be a good containment material for germanium. Confirmation under cyclic conditions 
has been demonstrated. 6 Germanium contained within graphite has proven to be a 
viable TES system. Furthermore, it has been demonstrated recently that this system 
virtually eliminates orbit-associated temperature variations in the storage systems. 7 The 
initial phase of this work has also demonstrated that germanium is compatible with 
pyrolytic carbon, CVD silicon carbide and CVD boron nitride. 

It has been demonstrated that NiSi-NiSi 2 can be readily produced and qualified. 
There appears to be a minor reaction between NiSi-NiSi 2 and pyrolytic carbon. However, 
the long-term compatibility of this system needs further investigation. CVD silicon 
carbide and CVD boron nitride also appear promising for containment of NiSi-NiSi 2 from 
chemical compatibility considerations. 
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